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A new measurement of the 6 — > Aj, production fraction multiplied by the A;, — > J/ip^ branching fraction was 
performed by the DO experiment using 6.1 fb~^ of pp collisions at y/s = 1.96 TeV. The result of this measurement, 
f{b At,)-B{At, J/tpJ^) = [6.01 ± 0.60 (stat.) ± 0.58 (syst.) ± 0.28 (PDG)] x IQ-^, represents an improvement in 
precision by about a factor of three with respect to the current world average. We give an estimate of B(At, J/ ip^)y 
which takes into account correlations among the different 6-hadron production fractions and other weakly decaying 
baryons. 



I — -i 1. Introduction 

Until recently, the only particle collider capable of producing b baryons was the Fermilab Tevatron Collider. Due 
to their relatively heavy mass, their production is suppressed with respect to the more favored B mesons, and 
even for the lightest and most copiously produced b baryon, the A{,{udb), only a few decay channels and properties 
have been studied. In particular, the uncertainties on Ab branching fractions are on the order of ^(30-60)%. With 
the full datasets of the DO and CDF experiments at Fermilab and the excellent performance of the Large Hadron 
Collider (LHC) and the experiments at CERN, it will be also possible to study more precisely important effects 
on b baryons, such as polarization, CP and T violation. 

The DO Collaboration reports a measurement of the production fraction multiplied by the branching fraction of 
the Ab — )• J/-0A decay relative to that of the decay i3° — J/^/jKg [1] , 

/(b ^ Aa) ■ g(Ab ^ J/M) ... 
" f{b ^ B") ■ B{BO ^ J/ijKl) ■ ^ ' 

Theestimationof /(&^ A,,)-S(A,, ^ J/V-A) is provided based on the best value of /(fa ^ 5°) -^(5" ^ J/i^Kl) [2]. 
A description of this analysis is given in the following sections. Finally, we give our estimate of S(Ab J/ipA). 



2. Experimental and theoretical status of B(Ab J/i/^A) 



The last measurement of (irci was performed by the CDF experiment [3] with only 7.8 ±3.4 A^ signal candidates. 
They found 



w.A. 

rcl 



0.27 ± 0.12 (stat.) ± 0.05 (syst.), (2) 



where W.A. stands for World Average. Based on this result the Particle Data Group (PDG) [2] reports 

f{b Ab) ■ B{Ab J/ipA) = (4.7 ± 2.3) x 10"^ (3) 



With more statistics and improved simulation of the processes involved and the experimental environment, this 
measurement can be greatly improved. It is important to mention that the dominant systematic uncertainty on 
this measurement is the unknown A^ polarization (section 5.3.1). 

On the other hand, there are several theoretical predictions of this branching fraction. For example, Ref. [4] uses 
perturbative QCD to find B{Ab J/V'A) - (1.65-5.27) x 10"''. The same branching fraction is calculated in the 
framework of the factorization hypotheses [5], using relativistic [6-8] and non-relativistic [9, 10] quark models, and 
ranges from ^ (1.1 — 6.1) x 10"'*. 
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a) 




Figure 1: Topology of the decays a) A(, J/tpK and b) i?" J/ipK%, with J/ip /i+p , A — > pvr and K% -k^-k 



The DO detector is described in detail in Ref. [11]. In general, in order to study B decays, the most relevant 
components are the central tracking system and the muon spectrometer. The DO central tracking system is 
composed of a silicon microstrip tracker (SMT) and a central fiber tracker (CFT) covering the pseudorapidity 
region \r}\ < 3.0 and \ri\ < 2.0, respectively, where 77 = — ln[tan(6'/2)] and 9 is the polar angle. They provide the 
ability to reconstruct charged tracks and vertices in a highly busy environment (typically more than 100 charged 
tracks) , and a surrounding 2 T superconducting solenoid allows precise measurements of the transverse momentum 
(pt) of the particles. The muon spectrometer consist of three layers of drift tubes and scintillator trigger counters, 
one located in front and two after 1.8 T iron toroids, and covering up to |ry| < 2.2. 



4. Data sample and event reconstruction 

This analysis uses an integrated luminosity of about 6.1 fb~^ recorded by the DO detector from 2002-2009 at 
-^i = 1.96 TeV. The data sample consists of events that satisfy single muon or dimuon triggers. 

4.1. Event reconstruction 

The decay topology of Ab — )■ J/ijjA and B'^ — )■ J/ipKg is shown in Fig. 1.^ The strategy to search for these 
decays is the following: 

(i) Look for events with two oppositely charged reconstructed muons, forming a common vertex, and with invari- 
ant mass Af(/i+/i^) in the range 2.8 — 3.35 GeV/c^. Muons are identified by matching tracks reconstructed 
in the central tracking system with track segments in the muon spectrometer. 

(ii) Search for pairs of oppositely charged tracks with a common vertex in those events satisfying the dimuon 
selection. For the A reconstruction, Monte Carlo (MC) studies support that the track with the highest pt 
is the proton. Events within 1.102 < M{pn-) < 1.130 GcV/c^ and 0.466 < M(7r+7r-) < 0.530 GeV/c^ are 



(iii) Ab and candidates are reconstructed by performing a constrained fit to a common vertex for the A or 
Kg candidate (a neutral track which is propagated from the pTr~ or tt'^'tt" common vertex according to the 
momentum direction of the pTT~ or tt+tt") and the two muon tracks. In this fit the dimuon mass is constrained 
to the W.A. J/V' mass [2]. Events within 5.0 < M{J/il;A) < 6.2 GeV/c^ and 4.8 < M{J/il)K%) < 5.8 GeV/c^ 
are selected. 



^Unless explicitly stated otherwise, the appearance of a specific charge state also implies its charge conjugate. 



3. Detector 



selected. 
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(iv) Finally, the reconstruction algorithm must be able to identify at least one pp interaction vertex^. In case 
of multiple interaction pp vertices in the event, the one closest to the B candidate vertex is tagged as the 
primary vertex (PV) for this candidate. 



4.2. Event quality and background suppression 

Several conditions are imposed on the quality of the reconstructed objects (tracks, vertices and parent particles): 

(v) Every muon track must be associated to at least two hits in (both) the SMT and CFT, and satisfy pT > 
2.0 GeV/c and \ri\ < 2.0. At least one muon must have segments in the muon system inside and outside the 
toroid. 

(vi) Since A and Kg are long-lived particles, they are likely to decay outside the beam pipe (and many of them 

outside the central tracking system). No detection condition is required in the SMT for the daughter tracks 
(proton or pion candidates); however, each of them must be detected with at least one hit in the CFT and, 
in total, they must not have more than two hits in the tracking detectors between the primary vertex and 
the common two-track vertex. Also, the impact parameter significance (the impact parameter with respect 
to the primary vertex divided by its uncertainty) is required to exceed 3 for both tracks and 4 for at least one 
of them. 

(vii) All (Ab, B^, A, Kg and J/tl>) decay vertices must be well reconstructed, with a probability greater than 
1%. 

In order to suppress undesirable backgrounds (distributed below or very close to the signal peaks) such as the 
cross- feed contamination^ between A;, and B^, and cascade decays of more massive baryons like E° A7 and 
An^, it is required that: 

(viii) Track pairs simultaneously identified as both A and Kg due to different mass assignments to the same tracks 
are removed. 

(ix) The pointing angle^ of the A (Kg) track to the J/ip vertex in the transverse plane must not exceed 2.5°. 

Finally, one can take advantage of the topology and kinematics of these decays in order to determine the final 
selection criteria. For example, it is easy to get rid of the prompt background (mainly J/i/j's coming from the 
primary vertex plus random tracks) by applying a minimum cut on the reconstructed decay length of the B 
particle. Similarly, the long-lived nature of the A and Kg can be used to suppress combinatorial background. 
To decide the final selection, MC events are generated for Af, — )• J/^p{ij,'^n~)A and B^ J/ip{ij,'^n~)Kg using 
PYTHIA [12] and evtgen [13] for the production and decay simulation, followed by full modeling of the detector 
response with GEANT [14], taking into account the effects of multiple interactions at high luminosity by overlaying 
hits from randomly triggered pp collisions on the digitized hits from MC, and event reconstruction as in data. Then 
the figure of merit S = Ng / \J Ns + Nb is maximized, where Ng is the number of signal candidates determined by 
MC and Nb the number of background candidates estimated by using data events in the sidebands of the expected 
signal. One ends up with the following requirements: 

(x) Dimuon candidates must satisfy pT{lJ.~^fi~) > 3.0 GeV/c. For the A {Kg), the pt must be greater than 1.6 
(1.0) GeV/c, the transverse decay length greater than 0.8 (0.4) cm and its significance greater than 4.0 (9.0). 
For the Af, (5°) candidate, the px must be greater than 5.0 GeV/c and the significance of the proper decay 
length^ is required to be greater than 2.0 (3.0). 



^The pp interaction vertex is determined by minimizing a funetion that depends on all reconstructed tracks in the event and a 
term that represents the average beam position constraint. 

^The Aj, sample may be contaminated with events that pass the A;, selection, or vice versa. 

*To be precise, this is the angle between the pT of the A and the vector from the J/V' vertex to the A decay vertex in the plane 
perpendicular to the beam direction. 

""The proper decay length is defined as LxyM/px, being px and M the transverse momentum and mass of the 6 hadron, respectively, 
and Lxy the distance between the primary vertex and the b hadron decay vertex in the transverse plane. 
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Figure 2: Mass distribution of a) Aft and b) B ^ candidates. The background distribution is parametrized as a second order 
polynomial and the signal distribution as a double Gaussian function. 



It may be the case that multiple candidates are found in the same event, for which only the one with the best 
probability of the B decay vertex is selected. 

The invariant mass distributions of the events satisfying the selection (i - x) are shown in Fig. 2. An unbinned 
likelihood fit to each distribution yields N/^^^j/^/^ = 314 ± 29 and Nga^j/^,xo = 2335 ± 73 candidates. 



5. Branching fraction measurement 
5.1. The method 

The number of observed A;, events (found in the previous section) decaying to J /ipK, with the J /tp going a pair 
of muons and A to a proton and a pion, is given by 

= Nprod [Ab ^ J/ij{n+ii-)K[pTT-)\ X eo [Ab ^ J ir)K{pT:-)\ , (4) 

where the number of decays produced in collisions is 

Nprod [Ab ^ J/^{^i+^l-)k{p^T-)\ = Ca{pp ^ hh) f{h ^ Ab)6 (Ab ^ J/^K) B {j/^ ^ n+ B {A ^ pn-) . (5) 

Here C is the integrated luminosity and a{pp — > bb) is the cross-section for the production of bb quarks. The 
detection efficiency eu [Ab J/tp{fi'^ fi^)A{pTT^)] encompasses acceptance effects as well as detector, trigger and 
reconstruction efficiencies for this decay. This efficiency is obtained from MC simulation. 

Similar expressions to Eqs. (4) and (5) can be obtained for — >■ J/'i{}{fi~^ ^~)Kg{n'^n^). Then, it is easy to 
show that 

^A.^J/^A g(gg ^ TT+TT-) 

'^ici — — • — TTT. r^; — • erci, (.oj 

Nbo^j/^k% B{A^ptt ) 

where dici is defined in Eq. (1) and the relative detection efficiency = £d[B'^ — !• J/^ifg]/e£)[Ab J/iPM is 
determined in the next section. For now, it is important to mention that most systematic and detector effects 
which are not fully implemented in the simulation (dead channels, trigger effects, pile-up, etc.) will cancel out in 
this ratio. Quantities such as b quark production, integrated luminosity and (to some extent) selection efficiencies 
are also canceled in ijrci. Hence the importance of choosing a normalization channel topologically equivalent to the 
decay under study. 



5.2. Detection efficiencies 

In order to determine the detection efficiencies, independent MC samples (different from the samples used to 
optimize the selection) of Ab and B^ decays are generated, with the same procedure described in section 4.2. 
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Proper Decay Length (K°) [cm] 



Figure 3: Proper decay length distributions for K% candidates reconstructed in the decay — >■ J/ipKs) comparing data 
and simulation. 



Important effects such as tracking detector efficiencies and luminosity dependence with time are incorporated in 
the simulation for different detector epochs and by the overlay of zero-bias events (triggered solely on the bunch 
crossing time). The same process for reconstructing and selecting events as in data is strictly followed. All the 
variables used in the selection are found to be in good agreement between data and MC. A noteworthy example is 
shown in Fig. 3, where the proper decay length distribution of Kg candidates is compared. Although no significant 
mismodeling was found in the simulation, any residual effect is expected to be reduced in the ratio of detection 
efficiencies. The relative detection efficiency of B'^ — > J/ipKg and Ab — > J/ipA decays is found to be 

erei = 2.37 ± 0.05 (MC stat.). (7) 

Using the number of reconstructed At, and B° decays in Fig. 2, B{K^ tt+tt") = 0.6920 ± 0.0005 and S(A 
pir-) = 0.639 ± 0.005 [2], we obtain 

cTrei = 0.345 ± 0.034 (stat.) ± 0.003 (PDG), (8) 
where the uncertainty due to inputs from the PDG [2] have been separated. 



5.3. Systematic uncertainties 

The sources of systematic uncertainties in the measurement of fXroi are summarized in Table I and explained 
below: 

1. The Af, and i3° yields can vary depending on the model used to describe signal and background in data and 
the mass range used in the analysis. No deviations larger that 5.5% with respect to the nominal value of CTicI 
are found in these tests. 

2. The cross- feed contamination between Af, and B'^ is quantified in simulation (see footnote 3). Due to this 
effect, the result on (Troi is estimated to change at most by 2.3%. 

3. The relative efficiency Crei depends on the models used in the simulation to decay the A;, and B^ particles. 

a) For B^, the SVSCP (scalar-vector-scalar with CP violation) method [13] in evtgen is used, resulting in 
a 2.0% deviation in a^c\- 

b) The Aft polarization can have a large effect on the A;, branching fraction. Since this is the dominant 
systematic uncertainty, we dedicate the following subsection to describe this phenomenon. 
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Table I: Systematics uncertainties on arc\- 



Source 


Error (%) 


Fit models 


5.5 


Cross-feed contamination 


2.3 


simulation 


2.0 


Ai, simulation (polarization) 


7.2 


Total (in quadrature) 


9.6 



5.3.1 . Ai, polarization 

Monte Carlo events with A;, initially polarized are generated following the methods used in [15]. The evtgen 
class HELAMP [13] was extended to accept one additional parameter that sets the value of the A^ polarization, P},. 
The polarization vector, given by 



P = 



z X p 



\z xp\' 

is set to the A(, particle through the spin density matrix 



Ph. 



p=l(l + a.p). 



(9) 



(10) 



The momentum p of the A^ particle is defined in the lab system, a are the Pauli Matrices and z is the direction 
of the proton beam. The HELAMP method decays Ab — )• J/'ip{ii'^ iJ.~)A{pTr~) according to the four complex helicity 
amplitudes, a± = A^±i^o ^^'^ ^± = -^tI.-I' ^^^^^ -Ma, A' denotes the amplitude for the A;, to decay into A and 

J/ip with helicities A and A'. The decay angular distribution depends on the angles 9 = {9,9i,62,4'i,4'2) depicted 
in Fig. 4a [16]. By integrating in four angles, it can be shown that only the 9 and <pi distributions depend on P;, 
{6 being the most relevant). In particular 9 follows the relation (see Fig. 4b), 



w{6 ;a±,b±, Pb) cx 1 + PbCtb cos 9, 
where the weak parity violating asymmetry parameter ai, is defined as 

|&+P-|a_|2_|6_|2 



■ 



K|2 + |?,+ |2 + |a_|2 + |6_| 



2 ■ 



(11) 



(12) 



Helicity amplitudes and polarization are independent unknown parameters which arc varied to study the effect 
on the A;, reconstruction efficiency. In particular, the slope of the cos^ distribution {abPb) is allowed to vary in 
the full range from -1 to 1. As expected, the largest variations are found in the extreme cases abPb = ±1, resulting 
in a 7.2% (maximum) deviation with respect to the nominal value of CTrei- 



6. Summary and discussion 

The relative production fraction times branching fraction of the A;, — >■ J/ipA to — >■ J/ipKg decays was 
measured using an integrated luminosity of 6.1 fb~^ collected with the DO experiment. The uncertainties in 
Eq. (8) can be combined in quadrature and the result, a^ei = 0.345 ± 0.047, can be compared with Eq. (2). 
The error is about 3 times smaller than in the previous measurement [3]. Equivalently, using the best value of 
f{b B°) ■ B{B° H> J/ipK^) = (1.74 ± 0.08) x 10"^ from the PDG [2], 



f{b Ab) ■ B{Ab J/V'A) = [6.01 ± 0.60 (stat.) ± 0.58 (syst.) ± 0.28 (PDG)] x 10"^ 

= (6.01 ±0.88) X 10-^ (13) 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



7 



a) 



A e 





P- 




» — 




' < 



0.1 
0.09 
0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




b) 



D0 Runll, MC 



MC Generation ( a^, Pj, = 1 ) 
MC Reconstruction 




P P 

Figure 4: a) At, — > J/i/;(/^"'"/i~)A(p7r~) decay angles and b) cos9 distribution of initially polarized At, with ObPt ~ 1. 



which can be compared with Eq. (3). 

The branching fraction B{Ab -> J /'^^) is slightly more difficult to report since there is not a published mea- 
surement of /(6 — >■ Kb). On the other hand, the DO and CDF experiments have observed other weakly de- 
caying baryons such as the S^, and il^, so the general assumption that /(6 — >■ btaryon) = fib — A;,) 
is not correct. A better approximation is to include the contribution of the S;, in the calculation, such that 
f{b — > bbaryon) ~ f [b — >■ Afo) -I" f {b Sj^) -|- /(& S°). Furthermore, we can assume isospin invariance to set 
f{b ■~-^,) = f{b Sg). It was also observed in Ref. [17] that f{b Bs)/f{b B°) « f{b ^ E^)/f{b Aft). 
Using the PDG values of f{b — )■ B°), f{b — ?► and f{b bbaryon) (from the combination of LEP and Tevatron 
results) and their correlations [2], we obtain 



B{Ab ^ J/VA) 



fjb^B^) 

fib b})(ij-yon) 



Jib^B, 



X BiB° JHjKI) X CTrel 



[11.08 ± 1.09 (stat.) ± 1.06 (syst.) ± 3.13 (PDG)] x 10"" 
(11.08 ±3.48) X lO"''. 



(14) 



The same assumptions on leads to B{h.b 



.w.A. K(K. ^ J/ipA) = (8.67 ±4.84) x lO"''. Both results are consistent 

within errors and favor theoretical models which predict a larger value for this branching ratio (see section 2). 

One final (but not less important) comment is that these measurements are useful to study & — >■ s decays such 
as Aft — > ^^fi^A. Due to their similar decay topology. Aft — > J/ipA can be used to normalize Aft /x+^^A: 



B(Aft ^ M+Ai"A) = 



..J /if! A 



X B{Ab -> J/V-A) X BiJ/ip ^ n+fi-) X 



eAb^j/i>A 



(15) 



The branching fraction of this rare decay, predicted to be ^ (2 — 5) x 10"® in the Standard Model [18-20], can be 
enhanced by new physics effects. Recent results by CDF report the observation of this decay [21]. Using the DO 
measurement of /(& Aft) • S(Aft J/ipA) [1], they found B{Ab ^J.~ A) = [1.73 ±0.42 (stat.) ±0.55 (syst.)] x 

10~® and no significant deviation from the Standard Model. 
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